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1 . 2  Work  in  a  180°  Curved  Passage  of  Circular  Cross-Section 

During  the  last  year  of  research  considerable  effort  was  expended 
on  obtaining  measurements  of  turbulent  flow  in  a  curved  pipe  and  its 
downstream  tangent.  The  main  objective  was  to  supply  Professor  Brian  Launder 
and  his  team  at  UMIST  with  data  useful  for  guiding  and  testing  numerical 
models  for  predicting  curved  pipe  flow. 

A  detailed  description  of  the  experimental  system,  the  measurement 
methodology,  the  results  and  their  discussion  Is  provided  in  the  Appendix. 
Listed  below  are  the  main  finding  of  the  study  and  the  recommendations 
for  future  work. 

Main  Findings 

1.  Laser-Doppler  measurements  of  the  longitudinal  and  circum¬ 
ferential  mean  velocity  components  and  their  normal  stresses  ware 
made  at  various  longitudinal  stations  in  a  180°  curved  pipe.  Special 
attention  was  given  to  documenting  the  flow  in  the  first  five 
hydraulic  diameters  of  the  downstream  tangent. 

2.  In  the  bend,  the  mean  longitudinal  velocity  component  changes 
little  after  a  90°  bend  angle,  but  the  circumferential  component 
never  achieves  a  fully-developed  state. 

3.  Large  levels  of  turbulent  stress  anisotropy  arise  everywhere 
In  the  bend  and  downstream  tangent. 

•  4.  Between  0  *  90°  In  the  bend  and  X/D  *  5  in  the  downstream 
tangent,  the  circumferential  velocity  profiles  display  a  reversal 
of  the  pressure-induced  secondary  flow. 

5.  Comparison  of  present  results  with  corresponding  measurements 
obtained  in  a  curved  duct  of  square  cross-section  (for  similar  flow 
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conditions)  shows  considerably  larger  levels  of  secondary  motion  in 
the  duct  of  square  cross-section. 

6.  Although  the  measurements  have  been  limited  to  a  specific 
4  -  plane,  they  should  serve  to  guide  and  test  theoretical 
developments  for  modeling  developing  turbulent  flow  in  curved  pipes. 
Recommendations 

1.  Further  experimentation  should  be  performed  to  measure  the 
extent  and  Intensity  of  the  extra  pair  of  counter-rotating  vortical 
structures  present  in  the  core  region  of  a  strongly  curved  pipe 

flow,.  In  particular,  the  relationship  between  these  structures  and  the 
local  levels  of  turbulence  Intensity  should  be  clarified. 

2,  Additional  and  more  extensive  measurements  are  required  to 
understand  the  flow  relaxation  process  taking  place  in  the  downstream 
tangent. 
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ABSTRACT 


Laser-Doppler  measurements  of  the  longitudinal  and  circumferential 
velocity  components  are  reported  for  developing  turbulent  flow  In  a 
strongly  curved  180°  pipe  and  Its  downstream  tangent.  In  the  bend,  the 
mean  longitudinal  velocity  component  changes  little  after  e  ■  90°,  but 
the  circumferential  component  never  achieves  a  fully-developed  state. 
Similar  behavior  Is  observed  In  the  normal  stresses,  with  large  levels  of 
flow  anisotropy  arising  everywhere  In  the  bend  and  downstream  tangent. 
Between  e  «  90°  and  X/D  ■  5,  the  circumferential  velocity  profiles 
display  reversals  of  the  secondary  flow  which  are  essentially  Independent 
of  the  Reynolds  number.  Ostensible  differences  between  turbulent  flows 
through  pipes  of  circular  and  square  cross-sections  are  revealed  by  com¬ 
paring  present  results  with  the  measurements  of  Chang,  et  al.  [5]. 
Although  limited  In  extent,  the  data  should  serve  to  guide  and  test 
theoretical  developments  for  modeling  developing  turbulent  flows  In 
curved  pipes. 


1.  INTRODUCTION 

1.1  The  Problem  of  Interest  and  Earlier  Work 

Over  the  last  ten  years,  approximately,  a  substantial  amount  of 
experimental  Information  has  been  obtained  for  developing  turbulent  flows 
In  curved  ducts  of  rectangular  cross-section  [1-7],  Detailed  measure¬ 
ments  In  these  Important  flow  configurations,  In  particular  of  the  secon¬ 
dary  motions  and  complex  turbulent  characteristics  which  they  Induce, 
have  been  made  possible  by  the  non-lntruslve  laser-Doppler  veloclmetry 
technique. 

In  contrast  to  the  rectangular  duct  configuration,  turbulent  flow- 
field  measurements  In  curved  pipes  appear  to  be  relatively  sparse  [8], 
This  Is  surprising  given  the  practical  Importance  of  the  curved  pipe  con¬ 
figuration.  The  lack  of  detailed  data  Is  partly  explained  by  the  greater 
difficulty  associated  with  optical  alignment  of  laser  beams  passing 
through  curved  surfaces  separating  media  of  different  Indices  of  refrac¬ 
tion. 

Unfortunately,  due  to  the  dominant  Influence  of  the  flat  walls  on 
the  secondary  motion  In  rectangular  ducts,  data  obtained  In  these  con¬ 
figurations  cannot  be  used  to  Interpret  accurately  the  characteristics  of 
flow  developing  In  curved  pipes.  For  example,  In  the  theoretical  analy¬ 
sis  by  Cuming  [9],  It  Is  shown  that  for  equal  radius  of  curvature  to 
diameter  ratio  and  dimensionless  streamwlse  pressure  drop,  the  ratio  of 
the  relative  Intensity  of  the  secondary  flow  In  a  duct  of  square  cross- 
section  to  that  In  a  pipe  of  circular  cross-section  Is  2.47  at  the  center 
line  location.  This  result  Is  strictly  valid  only  for  fully  developed 
laminar  flow.  Comparisons'  between  present  results  and  corresponding 


measurements  by  Chang,  et  al.  [5]  In  a  curved  duct  of  square  cross- 
section  with  similar  developing  flow  characteristics,  show  that  the 
secondary  motion  ratio  can  be  considerably  larger  than  2.47,  depending  on 
the  streamwlse  location.  This  point  will  be  quantified  more  precisely  in 
the  section  on  Results  and  Discussion. 

Among  the  early  measurements  of  secondary  turbulent  flow  patterns 
developing  In  and  downstream  of  curved  pipes  are  the  total  pressure  and 
yaw  results  obtained  by  Rowe  [10]  In  a  180°  bend.  With  reference  to  Fig. 
1,  the  longitudinal  stations  Investigated  were  e  ■  0®,  30®,  60®,  90°, 
120®,  150"  and  180®  In  the  bend,  and  X/D  ■  1,  5,  29  and  61  in  the 
downstream  tangent.  In  Rowe's  experiment  Rc/D  ■  12  and  Re  ■  2.36  x  10s, 
where  Rc  Is  the  bend  mean  radius  of  curvature,  D  its  diameter  and  Re  Is 
the  Reynolds  number  based  on  the  bulk  fluid  velocity.  An  upstream 
tangent  of  length  X/D  ■  69  was  used  to  provide  a  fully  developed  axl sym¬ 
metric  turbulent  velocity  profile  at  the  bend  Inlet  plane.  Among  Rowe's 
findings  were  that:  1)  secondary  flow  Is  most  Intense  at  about  d  ■  30®, 
as  of  which  point  the  total  pressure  gradient  Induces  a  streamwlse  com¬ 
ponent  of  vortlclty  opposite  In  sense  of  rotation  to  the  streamwlse  vor- 
tlclty  produced  at  the  start  of  the  bend;  2)  the  curved  flow  Is 
essentially  fully  developed  past  0  *  90°;  3)  In  the  downstream  tangent 
there  Is  a  gradual  reversion  to  fully  developed  turbulent  straight  pipe 
flow  with  the  secondary  flow  persisting  past  X/D  •  29  but  being  undetec¬ 
table  at  X/D  -  61;  4)  there  Is  evidence  of  local  reversal  In  the  secon¬ 
dary  flow  direction  along  the  bend  symmetry  plane  between  e  *  90®  and 
X/D  “  5;  5)  the  secondary  flow  appears  to  be  most  concentrated  at  the 
bend  Inner  radius  of  curvature;  6)  at  all  longitudinal  stations,  the 
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maximum  In  the  streamwlse  velocity  component  was  always  located  between 
the  bend  center  line  and  the  outer  radius  wall. 

Because  of  the  nature  of  the  variables  measured  In  Rowe's  study,  and 

because  of  the  Influence  of  the  mechanical  probes  used  In  the  flow,  an 

*• 

accurate  determination  of  the  fully  developed  flow  streamwlse  location 
was  not  possible.  A  combination  of  weak  secondary  motions,  undetected  by 
the  yawmeter,  with  unchanging  relatively  large  total  pressure  measure¬ 
ments  would  certainly  have  given  the  Impression  of  fully  developed  flow. 
Due  to  experimental  uncertainty,  Rowe  was  not  able  to  define  precisely 
the  magnitude  and  extent  of  secondary  flow  reversal  In  the  vicinity  of 
the  bend  symmetry  plane.  Although  valuable  for  helping  to  understand  the 
behavior  of  the  mean  flow,  particularly  In  relation  to  Invlscld  effects, 
the  pressure  and  yaw  contours  obtained  by  Rowe  shed  no  light  on  the  tur¬ 
bulent  characteristics  of  the  flow. 

In  the  only  laser-veloclmetry  study  known  to  us  for  developing  tur¬ 
bulent  flow  In  curved  pipes,  Enayet,  et  al.  [11]  Investigated  the  motion 
In  a  90°  bend  with  Rc/D  »  2.8  and  Re  ■  43,000.  In  their  study  relatively 
thin  boundary  layers  were  Induced  at  the  bend  Inlet  plane  by  means  of  a 
smooth  contraction.  The  measurements  consisted  of  longitudinal 
(streamwlse)  components  of  mean  and  fluctuating  velocities  only,  supple¬ 
mented  by  wall -pressure  measurements.  As  expected,  the  measurements 
suggest  the  development  of  a  strong  secondary  flow  In  thr>  form  of  a  pair 
of  counter-rotating  longitudinal  vortices.  The  secondary  flow  charac¬ 
teristics  were  found  to  depend  on  the  thickness  of  the  Inlet  boundary 
layers.  In  addition,  the  relatively  flat  Inlet  velocity  profile  and 
small  value  of  Rc/D  in  their  experiment  enhanced  streamwlse  acceleration 
of  the  flow  at  the  Inner  radius  wall  relative  to  the  outer  wall.  As  a 


result,  throughout  the  bend  and  up  to  X/D  ■  6  in  the  downstream  tangent, 
the  peak  in  the  streamwise  velocity  component  remained  between  the  bend 
center  line  and  the  Inner  radius  wall.  This  finding  contrasts  strc  y 
with  finding  (6)  of  Rowe's  study.  Because  the  Dean  numbers  of  these  two 
flows  are  very  similar,  being  De  ■  Re  (D/Rc)1^2  «  6.8  x  10^  In  Rowe's 
bend  and  2.6  x  10^  In  the  bend  of  Enayet,  et  al.,  the  difference  In  the 
maximum  velocity  location  must  be  due  to  the  different  Inlet  flow  boun¬ 
dary  layer  thicknesses.  A  similar  finding  has  been  reported  In  [3]  for 
the  flow  through  a  curved  duct  of  square  cross  section. 

In  addition  to  the  data  taken  at  three  longitudinal  stations  In  the 
curved  pipe,  (a  »  30®,  60®  and  90®),  Enayet,  et  al.  [11]  made  measure¬ 
ments  at  X/D  ■  1  and  6  In  the  downstream  tangent.  Although  the  authors 
did  not  measure  the  cross-stream  velocity  components  at  any  of  these  sta¬ 
tions,  they  reported:  "a  large  secondary  flow  ...  now  extends  over  the 
entire  flow  area  and  persists  In  the  plane  six  diameters  downstream  of 
the  band  ..."  In  the  absence  of  secondary  component  measurements,  the 
Inference  concerning  secondary  flow  In  the  downstream  tangent  could  only 
have  been  based  on  the  distorted  appearance  of  the  streamwise  velocity 
component  contours  at  X/D  «  1  and  6,  respectively. 

1.2  The  Present  Contribution 

The  above,  admittedly  brief,  review  points  to  a  serious  deficiency 
In  the  availability  of  accurate  experimental  information  for  Illuminating 
the  physics  of  developing  turbulent  flow  through  curved  pipes.  It  Is 
also  clear  that  little  seems  to  be  known  about  the  relaxation  process 
undergone  by  fluid  leaving  a  curved  pipe  to  enter  a  straight  pipe  sec¬ 
tion.  The  situation  Is  especially  serious  since  for  this  class  of 
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complex  three-dimensional  flows  such  data  Is  Indispensable  for  guiding 
the  development  and  testing  of  theoretical  models  and  calculation  proce¬ 
dures. 

The  purpose  of  this  work  has  been  to  obtain  limited  but  critical 
measurements  of  the  longitudinal  and  circumferential  mean  velocity  com¬ 
ponents  and  their  respective  turbulence  Intensities  at  selected  longitu¬ 
dinal  stations  in  a  180®  curved  pipe  and  Its  downstream  tangent.  In 
addition  to  Its  academic  significance,  this  configuration  Is  especially 
relevant  to  heat  exchangers  and  particle-laden  flow  equipment.  In  the 
study,  special  emphasis  has  been  placed  on  documenting  the  flow  between 
e  ■  177°  In  the  bend  and  X/D  ■  5  In  the  downstream  tangent. 

2.  THE  EXPERIMENT 

2.1  Flew  System  and  Instrumentation 

The  experimental  system  was  composed  of:  a  water  rig,  of  which  the 
most  Important  component  was  the  flow  test  section;  a  laser-Doppler  velo- 
clmeter  and  Its  associated  electronic  Instrumentation;  and  a  POP  11/34 
Digital  Equipment  Corporation  minicomputer. 

•  The  basic  components  of  the  flow  test  section  are  shown  schemati¬ 
cally  In  Fig.  1,  They  comprised  two  straight  pipes  and  a  180°  curved 
pipe,  constructed  from  transparent  plexiglass.  The  pipe  cross-section 
was  circular  throughout  with  a  4.45  cm  ±  0.02  cm  Inner  diameter  (0).  The 
ratio  of  bend  mean  radius  of  curvature  (Rc)  to  pipe  diameter  was 
Rc/D  ■  3.375,  Both  tangents  were  of  length  x  ■  54.7  D,  being  respec¬ 
tively  attached  to  the  0°  (Inlet)  and  180°  (outlet)  planes  of  the  bend  by 
means  of  flanges.  Special  care  was  taken  to  avoid  possible  mismatches 
between  the  component  cross-sections  which  otherwise  might  disturb  the 
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flow.  A  stainless  steel  (20  mesh)  screen  and  a  3.5  cm  long  flow 
straightening  honeycomb  section  were  placed  upstream  of  the  entrance 
tangent.  They  were  held  In  place  by  a  thin  plexiglass  plate.  The  plate 
was  3.175  mm  thick  and  had  85  holes  of  3.175  mm  diameter  arranged  in  a 
rectangular  array  spaced  4.495  mm  on  the  centers  In  each  direction.  The 
purpose  of  the  screen-honeycomb-plate  combination  was  to  make  uniform  the 
cross-stream  plane  distribution  of  the  flow  an,  .ccelerate  Its  streamwlse 
development. 

The  curved  pipe  section  was  constructed  by  fitting  together  two  sym¬ 
metrical  half  sections  of  plexiglass,  each  respectively  machined  on  one 
of  Its  flat  faces  to  contain  the  shape  of  a  semicircular  open  channel. 
This  method  of  construction  ensured  that  when  matched  at  the  common  sym¬ 
metry  plane  the  cross-section  of  the  resulting  curved  pipe  was  accurately 
circular.  When  assembled,  the  curved  pipe  section  had  flat  outer  sur¬ 
faces  of  thickness  d  ■  1.48  ±  0.02  cm  (see  Fig.  1).  The  straight  pipe 
sections  had  a  circumferentially  uniform  wall  thickness  of  0,32  ±  0.02 
cm.  The  entire  test  section  was  supported  horizontally  by  a  uni  struct 
structure,  with  the  bend  symmetry  plane  parallel  to  the  floor. 

The  test  section  was  part  of  a  closed  loop  through  which  water  was 
made  to  flow  by  gravity  from  a  constant  head  tank.  From  this  tank  the 
flow  passed  through  the  test  section,  and  then  Into  a  large  sump  tank 
from  where  It  was  pumped  back  to  the  constant  head  tank.  Flow  to  the 
head  tank  was  controlled  by  a  gate  valve  and  measured  using  a  venturi 
meter  connected  to  a  differential  mercury  manometer.  The  possibility  of 
propagating  disturbances  through  the  test  section  was' eliminated  by 
avoiding  the  use  of  valves,  sharp  bends  and  metering  devices  In  the  test 
section  flow  loop.  Flow  to  and  from  the  test  section  tangents  was  chan- 
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neled  through  2  In  (5,08  cm)  l.d.  gently  bent  tygon  tube  pieces.  Baffles 
and  screens  located  In  the  constant  head  tank  served  to  dampen  the 
swirling  motion  of  the  flow  leaving  the  tank.  Residual  swirl  and  weak 
secondary  motlbns  Induced  by  the  mild  curvature  In  the  tygon  tube 
upstream  of  the  test  section  were  eliminated  by  the  flow  straightening 
section. 

Measurements  of  the  mean  flow  and  turbulence  characteristics  were 
made  using  a  DISA  55X  modular  series  laser-Doppler  veloclmeter  In 
backscatter  mode.  The  system  has  already  been  described  In  full  detail 
In  reference  [5].  A  2-watt  Lexel  Argon-Ion  water-cooled  laser  was  used 
as  the  light  source.  The  laser  and  veloclmeter  optics  were  mounted  to 
the  top  of  a  thick  aluminum  table  which  was  Itself  firmly  bolted  to  an  x, 
y,  z  traversing  mechanism.  The  traversing  mechanism  cou.ld  displace  the 
table  top  *  7.5  cm  In  5  um  Increments  along  any  of  the  coordinate  axes  by 
means  of  three  linearly  encoded  stepping  motors  monitored  by  the  POP 
11/34  minicomputer.  The  minicomputer  functioned  as  the  central  data 
acquisition  and  reduction  controller  of  data  validated  and  measured  by  a 
DISA  55L96  Doppler  signal  procesor  or  "counter."  Directional  ambiguity 
In  the  circumferential  (or  tangential)  velocity  component  was  resolved  by 
means  of  a  Bragg  cell  combined  with  electronic  downmixing. 

2.2  Experimental  Methodology 

Prior  to  an  experimental  run,  water  was  allowed  to  flow  through  the 
rig  until  It  was  purged  of  air  bubbles  and  had  attained  a  steady  thermal 
state  corresponding  to  about  20°C.  Mass  flow  through  the  test  section 
was  controlled  by  setting  the  constant  head  overflow  condition  to  a  mere 
trickle  and  continuously  monitoring  the  pressure  drop  through  the  venturi 
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meter  connected  to  the  head  tank  feed  line.  Experiments  were  performed 
for  two  values  of  the  bulk  average  velocity  corresponding  to 
Ujj  »  1.29  ±  0.03  m/s  and  ■  2.47  ±  0.04  m/s,  respectively.  During  the 
course  of  an  experimental  run  the  fluid  temperature  rose  from  18°C  to 
22°C,  approximately.  Values  for  the  Reynolds  numbers  of  these  two  flows, 
based  on  physical  properties  of  water  at  208C,  are  Re  -  57,400  and 
110,000,  respectively.  The  associated  Dean  numbers  are  De  •  31,300  and 
59,900. 

The  veloclmeter  optical  probe  volume  was  formed  by  the  Intersection 
of  two  514.5  nm  (green)  light  beams  with  a  half-angle  In  air  of  3. 28°. 

The  measurement  volume  charaterl sties  were:  a  diameter  of  0.1  mm,  a 
length  of  1.6  nm,  and  a  fringe  spacing  of  4.50  urn  with  about  18  fringes 
contained  In  the  optical  probe.  In  reality  spatial  filtering  and 
threshold  settings  on  the  counter  reduced  the  dimensions  of  the  optical 
probe.  The  probe  volume  was  positioned  by  fine  control  of  the  motorized 
traversing  table. 

At  each  measurement  station  Initial  positioning  of  the  laser  beams 
with  respect  to  the  plexiglass  surface  was  achieved  through  the  use  of 
*  reference  markers.  Paper  drawn  with  precise  cross  hairs  was  taped  onto 
the  flat  outer  surface  of  the  curved  pipe  section  and  the  beams  were 
aligned  with  respect  to  the  marks,  checking  that  no  sideways  beam 
deviation  occurred  as  the  beams  were  displaced  In  the  vertical  direction. 
After  the  paper  was  removed,  symmetrical  reflections  facilitated  vertical 
alignment  of  the  optical  volume  In  the  cross-stream  plane.  Alignment  of 
the  beams  In  the  straight  pipe  sections  was  performed  In  a  similar 
manner.  Special  care  was  taken  to  ensure  that  the  optical  probe  fringe 
pattern  was  always  aligned  normal  to  the  velocity  component  of  Interest 
to  within  ±  0.3°. 


With  a  reference  position  established,  usually  on  the  outer  surface 
of  the  test  section,  the  computer  software  was  activated  which  controls 
signal  acquisition  and  data  processing.  An  Initial  check  on  flow  sym¬ 
metry  was  performed  at  X/D  ■  -2  In  the  upstream  tangent.  Measurements  of 
the  streamwlse  velocity  component  and  its  normal  stress  at  this  station 
(half  profiles  are  shown  In  Figs.  2  and  4)  are  In  good  agreement  with  the 
data  of  Laufer  [19].  Measurements  of  the  circumferential  velocity  com¬ 
ponent  (not  plotted)  showed  that  It  was  everywhere  less  than  -  1.5%  of 
the  bulk:  average  velocity  at  this  station.  Calculations  of  U^,  obtained 
by  Integrating  the  measured  velocity  distributions  at  X/D  ■  -2,  yielded 
values  In  good  agreement  (to  within  -  3%)  with  the  values  obtained  from 
the  venturi  meter  pressure  drop  readings. 

Following  the  symmetry  and  mass  flow  confirmations,  all  subsequent 
measurements  wore  restricted  to  vertical  (radial)  scans  In  the  plane 
+  «*/2  In  the  symmetrical  upper  half  of  the  test  section.  Scans  were 
made  at  X/D  ■  -2,  -1,  1,  2,  3,  4  and  5  In  the  straight  pipes  and  e  ■  30°, 
45°,  90°,  135°,  and  177°  In  the  bend.  At  each  of  these  stations  15 
radial  positions  were  probed,  starting  at  the  centerline  and  moving  In 
increments  of  1.5  mm  toward  the  wall,  with  a  16th  position  fixed  at  0,73 
mm  from  the  Inside  pipe  wall.  At  each  measurement  point  the  mean  flow 
and  turbulence  characteristics  were  statistically  determined  from  popula¬ 
tions  of  1  to  3  samples  consisting  of  1,000  Individual  measurements  each. 
More  than  90%  of  all  the  measurements  consisted  of  at  least  2  samples  per 
point.  Each  Individual  measurement  was  required  to  satisfy  the  counter 
5/8  validation  comparison  to  within  a  preset  tolerance  of  6%.  At  every 
validation  of  a  Doppler  burst  a  "data  ready"  signal  was  Issued  by  the 
counter  to  a  logic  conversion  circuit.  This  circuit  then  sent  a 


triggering  pulse  to  the  computer  parallel  line  Interface  module  which  was 
checked  for  data  availability  by  a  software  loop  approximately  every  20 
vs. 

Attempts  were  made  to  optimize  signal  quality  and  particle  con¬ 
centration.  Corn  starch,  with  particles  ranging*  In  size  between  1  and  10 
um*  was  found  to  give  a  consistently  high  slgnal-to-nolse  ratio  In  addi¬ 
tion  to  a  low  settling  rate.  Alternative  particle  types  rejected  Include 
T102,  non-dairy  creamers,  and  mother-of-pearl.  Data  rates  varied  from 
50-1000  Hz,  with  counter  validation  levels  ranging  between  10  and  90% 
approximately.  Factors  Influencing  these  levels  Included,  in  order  of 
relative  Importance:  optical  alignment,  thickness  of  the  wall,  probe 
volume  location  In  the  flow  field,  and  seeding  concentration. 

Although  the  veloclmeter  Is  capable  of  two-component  measurements, 
curvature  of  the  test  section  surfaces  precluded  dual-channel  operation. 
Values  for  the  streamwlse  component  of  tne  mean  and  fluctuating 
velocities,  U  and  u',  were  derived  from  measurements  obtained  with  the 
probe  volume  fringes  aligned  normal  to  the  longitudinal  coordinate  direc¬ 
tion.  Values  for  the  circumferential  components,  V  and  v',  were  obtained 
through  a  rotation  of  the  optics  by  90°,  thus  placing  the  fringes  normal 
to  the  circumferential  component  direction. 

Because  of  the  large  positive  magnitudes  of  the  streamwlse  velocity 
component,  frequency  shifting  was  unnecessary  for  measurements  of  this 
component.  By  contrast,  the  circumferential  velocity  component  exhibited 
low  magnitudes,  with  signs  both  positive  and  negative.  To  resolve  this 
directional  ambiguity,  a  net  frequency  shift  of  500  kHz  was  Imposed  using 
a  0ISA  55N10  Bragg  cell  combined  with  electronic  downmixing.  Appendix  II 
provides  an  explanation  for  Implementing  the  frequency  shift  procedure. 
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2.3  Measurement  Error 

Error  sources  affecting  the  accuracy  (systematic  error)  and  preci¬ 
sion  (random  error)  of  laser-Doppler  measurements  have  been  discussed  by, 
for  example.  Durst,  et  al.  [12],  Drain  [13]  and  Buchhave  [14].  In  this 
study  the  most  serious  potential  systematic  errors  Mere  attributed  to 
probe  volume  positioning,  velocity  gradient  broadening  and  velocity  bias, 
respectively.  The  main  sources  of  random  error  affecting  the  precision 
of  the  measurements  were  due  to:  a)  statistical  sampling  uncertainty  (due 
to  the  finite  size  of  sample  populations);  and  b)  uncertainty  In  the 
determination  of  the  normalizing  bulk  average  velocity,  Ufa.  Estimates  of 
the  first  uncertainty  ware  derived  from  the  measurements  themselves.  For 
both  mean  velocity  components,  this  uncertainty  was  always  less  than 
*  2-3X  (r.m.s.  error),  while  for  the  turbulence  Intensities  It  was  always 
less  than  t  2-5%.  The  error  In  was  *  2%  and  arose  principally  from 
uncertainties  In  the  manometer  pressure  readings. 

The  Initial  probe  volume  location  uncertainty  was  estimated  to  be 
half  the  effective  probe  volume  length,  or ±  0.5  mm  approximately.  This 
error  Is  altered  as  the  beams  traverse  the  pipe  because  of  Index  of 
refraction  and  wall  curvature  effects.  To  calculate  the  actual  probe 
volume  location,  and  to  estimate  the  uncertainty  In  this  location,  It  Is 
necessary  to  know  the  relevant  trigonometric  relations.  Two  of  the  three 
cases  of  Interest  here  have  been  analyzed  by  Blcen  [15]  and  Boadway  and 
Karahan  [16],  respectively.  The  present  authors  have  found  Inconsisten¬ 
cies  In  the  analyses  In  these  references,  as  well  as  between  them  for 
Identical  configurations.  For  this  reason  It  was  conslcered  prudent  to 
rederive  the  trigonometric  relations  of  relevance  to  this  study. 

Appendix  I  presents  relations  corresponding  to  flows  bounded  by:  a) 
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parallel  flat  surfaces;  b)  concentric  cylindrical  surfaces;  and,  e)  one 
flat  and  one  cylindrical  surface. 

For  parallel  and  concentric  cylindrical  surfaces,  references  [15-16] 
obtained: 

U  .  CfUM 

In  this  relation  Is  the  velocity  measurement  obtained  In,  for  example, 
water  using  a  half  angle  between  beams  In,  for  example,  air.  The  factor 
Cf  is  necessary  to  “correct"  the  measurement  so  that  the  result  repre¬ 
sents  the  water  velocity,  for  parallel  surfaces,  to  a  small  angle 
approximation  It  can  be  shown  that: 

U  ■  (a/r )  ^  "  Hf  Um 

where:  a  Is  the  beam  half  angle  In  air,  y  Is  the  beam  half  angle  In 
water,  andnf  Is  the  Index  of  refraction  of  water  relative  to  air  (here 
taken  as  1.333).  Curiously,  reference  [15]  takes  Cf  ■  1  for  parallel 
surfaces,  a  result  valid  only  for  air  as  the  working  fluid.  Similarly, 
the  expression  In  [16]  for  for  correcting  circumferential  velocity 
measurements  Is  In  error  by  a  factor  of  r^. 

Table  1  summarizes  the  relations  derived  for  the  various  cases  ana¬ 
lyzed  In  Appendix  I.  The  relations  depend  on  the  velocity  component  of 
Interest  and  the  shape  of  the  surfaces  at  the  measurement  station  In 
question.  For  example,  the  beam  orientation  required  to  measure  the 
longitudinal  velocity  component  In  the  plane  $  -ir/2  renders  the  beams 
Insensitive  to  surface  curvature  effects.  For  this  component  the  beams 
are  crossed  such  that  their  bisector  Is  normal  to  the  test  section  sur¬ 
face  and  the  beams  are  contained  In  the  4>  *ir/2  plane.  Therefore,  for 
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all  practical  purposes,  both  In  the  curved  and  straight  pipe  sections  the 
beams  behave  as  If  they  were  passing  through  a  pair  of  parallel  Inter¬ 
faces^.  By  contrast  the  orientation  required  to  measure  the  circumferen¬ 
tial  velocity  component  makes  It  necessary  to  distinguish  between  the 
straight  and  curved  pipe  sections.  For  this  component  the  beams  are 
crossed  such  that  their  bisector  Is  normal  to  the  test  section  surface 
and  Is  contained  In  the  $  “ir/2  plane,  while  the  beams  themselves  are 
contained  In  a  plane  perpendicular  to  the  +  «*/2  plane. 

Table  1  shows  that  In  regions  of  surface  curvature,  both  the  probe 
volume  location  and  the  half  angle  between  the  beams  In  water  depend  non- 
linearly  on  the  probe  volume  virtual  radial  location  (the  location  In 
air).  Calculations  of  the  error  Incurred  Iny.  the  half  angle  of  the 
beams  In  water,  due  to  at  0.5  mm  positioning  error  In  the  Initial  probe 
volume  location  showed  that  this  source  of  Inaccuracy  was  negligible. 
Similarly,  calculations  showed  that  the  Initial  positioning  uncertainty 
was  only  significant  very  near  the  pipe  walls,  In  regions  with  steep  gra¬ 
dients.  Thus,  estimates  of  the  maximum  (absolute)  systematic  errors  In 
the  measurements  obtained  nearest  the  walls  gave:  -  4%  for  U/Ub;  -  10% 
for  u'/Ub;  -  15%  for  V/Ub,  and-.  10%  for  v ' /Ub*  In  the  bulk  of  the  flow 
all  these  errors  were  significantly  less  than  -  2%. 

Velocity  gradient  broadening  has  been  analyzed  by  Melling  [18]  who 
proposed  a  simple  method  for  estimating  Its  magnitude.  For  the  con- 


******************** 

^Due  to  the  small  angle  between  the  beams,  the  error  Incurred  in  the 
bend  as  a  consequence  of  making  this  approximation  Is  negligible. 
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dltlons  of  this  study,-  gradient  broadening  In  the  measurements  was  esti¬ 
mated  to  be  less  than  1%  everywhere  In  the  flow. 

Various  weighting  methods  have  been  proposed  to  correct  for  the 
velocity  bias  effect  [13,14],  but  none  of  these  Is  entirely  satisfactory. 
They  all  Involve  assumptions  regarding  the  statistical  distribution  of 
particles  In  the  flow.  In  addition,  the  corrections  can  be  Influenced  by 
the  problem  of  ".Incomplete-signal  bias."  The  use  of  the  weighting  scheme 
proposed  by  McLaughlin  and  Tlederman  [17],  In  which  the  particle  resi¬ 
dence  time  Is  taken  to  be  Inversely  proportional  to  the  modulus  of  the 
measured  velocity,  was  investigated  In  this  study.  For  the  longitudinal 
velocity  component,  differences  between  weighted  and  unweighted  oata  sets 
were  relatively  small  (Jess  than  1%).  For  the  circumferential  component 
the  differences  were  much  larger,  with  the  weighted  data  giving  a  less 
realistic  (overdamped)  picture  of  the  flow  field.  The  reason  for  this  Is 
attributed  to  an  overcorrection  In  the  weighting  scheme  which,  for  the 
circumferential  component,  amounts  to  a  disproportionate  weighting  of 

i 

slow  moving  particles.  This  occurs  as  a  consequence  of  taking  the  par¬ 
ticle  residence  time  as  tres  -  V"1  Instead  of  tpes  *  (U2  +  V2  +  W2)"1/r2. 
The  latter,  more  accurate,  weighting  could  not  be  performed  with  the  pre¬ 
sent  single  component  measurement  system. 

For  turbulence  Intensities  less  than  10%  Drain  [13]  shows  that  velo¬ 
city  bias  corrections  of  the  mean  and  r.m.s.  velocities  should  not  amount 
to  more  than  -  2%.  In  this  study  the  measured  turbulence  Intensities 
never  exceed  -  12%,  As  a  result,  plots  of  the  unweighted  data  only  are 
presented  here.  Appendix  III  contains  a  listing  of  both  the  weighted  and 
unweighted  measurements. 
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3.  RESULTS  AND  DISCUSSION 

Measurements  of  the  longitudinal  (U/Ub)  and  circumferential  (V/Ujj) 
velocity  components  and  their  respective  turbulence  Intensities  (u'/Ub, 
v'/Ub)  are  shown  In  Figs.  2*5  and  are  tabulated  In  Appendix  III.  Except 
for  relatively  small  differences,  the  dimensionless  profiles  show  that 
flow  development  Is  essentially  Independent  of  Reynolds  number  for 
Re  >  57,400. 

3.1  The  Inlet  Flow 

The  U  and  u'  component  measurements  at  X/D  *  -2  show  good  agreement 
with  the  developed  flow  measurements  obtained  by  Laufer  [19]  In  a 
straight  pipe.  At  X/D  •  -2  the  V  component  Is  less  than  (+)  0.015 
Ub  (not  plotted).  However,  at  X/D  -1  and  6  *  3°,  the  V  component  Is 
-  {+)  0,04  Ub.  with  reference  to  the  coordinate  system  in  Fig.  1,  this 
implies  that  the  Inlet  flow  Is  gradually  (but  only  slightly)  accelerated 
towards  the  Inner  radius  of  the  pipe,  a  result  expected  from  potential 
flow  theory  and  previously  observed  by  others. 

3.2  Mean  Flow  and  Turbulence  Intensity 

In  the  first  half  of  the  bend  (e  «  3®  to  90®),  the  U-proflles  show 
the  core  of  the  streamwise  flow  losing  speed  while  the  flow  near  the  wall 
accelerates.  V-profllas  at  corresponding  stations  reveal  the  development 
of  a  strong  secondary  flow.  This  secondary  flow  is  induced  by  the  trans¬ 
verse  pressure  gradient  set  up  between  the  outer  (rQ)  and  inner  (r^ )  wall 
regions  of  the  bend.  In  the  pipe  center,  it  works  to  overcome  and 
reverse  the  sense  of  the  cross-stream  motion  in  the  Inlet  flow.  Near  the 


wall  the  original  sense  of  cross-stream  motion  Is  preserved  and  Inten¬ 
sified,  with  (V/Ub)max  -  0.30  at  9  -  45®. 

Between  a  ■  45®  and  135°  the  V-proflles  reveal  a  striking  feature  of 
the  flow.  In  the  fluid  core,  r/(0/2)  <  0.5,  the  cross-stream  flow 
undergoes  a  second  reversal  In  its  sense  of  motion  and  Is  redirected 

a 

towards  the  pipe  Inner  radius.  As  a  result,  a  region  of  negative  cross¬ 
stream  flow  (directed  from  r^  towards  rQ)  Is  trapped  between  the  core  and 
the  wall;  for  example,  V/Ub  ~  -0.05  at  e  ■  135*.  Between  e  ■  135°  In  the 
bend  and  X/D  *  5  In  the  downstream  tangent,  the  region  of  negative  cross¬ 
stream  flow  gradually  disappears.  Notwithstanding,  an  Imprint  of  Its 
presence  remains  In  all  the  V-proflles. 

The  occurrence  of  a  second  cross-stream  flow  reversal  past  e  «  90® 

In  the  bend  supports  the  concept  proposed  by  Rowe  [10],  that  each  sym¬ 
metrical  half  section  of  the  bend  develops  two  counter-rotating  vortical 
structures.  The  more  Intense  of  the  two,  located  between  the  pipe  wall 
and  the  core  of  the  flow,  preserves  the  sense  of  cross-stream  motion 
Induced  by  the  transverse  pressure  gradient  at  the  start  of  the  bend. 

The  smaller,  weaker,  structure  Is  mainly  confined  to  the  core  and  Is 
attributed  to  the  formation  of  a  transverse  pressure  gradient  opposite  In 
sign  to  that  at  the  start  of  the  bend. 

The  present  mean  flow  measurements  place  on  a  firm  quantitative 
basis  the  Interpretation  by  Rowe  of  his  total  pressure  and  yaw  measure¬ 
ments.  However,  In  contrast  to  Rowe's  results,  the  V-proflles  provided 
here  show  that  fully  developed  mean  flow  conditions  are  not  established 
In  the  bend. 

There  Is  a  qualitative  similarity  between  the  present  V-proflles  and 
corresponding  laminar  flow  measurements  made  by  Olson  and  Snyder  [20]  In 
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300°  curved  pipes.  In  their  experiments  two  curvature  ratios 
(Rc/D  *  2.33  and  8)  were  Investigated  for  two  values  of  the  Reynolds 
number  (Re  ■  300  and  1030).  The  Inlet  flow  was  contracted  to  produce 
nearly  flat  entrance  profiles.  Measurements  of  the  cross-stream  flow 
were  made  at  several  longitudinal  stations  with  'a  pulsed-probe  anemo¬ 
meter.  Between  e  ■  40°  and  90°  the  authors  found  a  cross-stream  flow 
reversal  In  the  center  of  the  pipe  similar  to  the  present  one. 

The  longitudinal  variation  of  the  local  ratio  V/U  Is  plotted  In  Fig. 
6  for  two  values  of  r/(D/2)  in  the  present  flow.  The  profiles  show 
clearly  the  cross-stream  flow  reversals  which  arise  in  the  bend.  They 
also  show  that,  although  substantially  reduced,  secondary  motion  persists 
between  e  ■  177°  and  X/D  «  5  In  the  downstream  tangent,  and  Is  stronger 
near  the  pipe  center  than  at  the  wall. 

Figure  6  also  summarizes  the  most  Interesting  variations  of  the  tur¬ 
bulence  components  u'/Ub  and  v'/Ub  in  the  bend  and  downstream  tangent. 
Actually  plotted  Is  the  ratio  v'/u',  a  measure  of  turbulent  flow  ani¬ 
sotropy.  Because  the  data  are  relatively  few  and  restricted  to  the 
4  *  ir/2  plane  shown  In  Fig.  1,  and  because  the  physical  processes 
governing  the  generation,  redistribution  and  destruction  of  the  turbulent 
stresses  are  complex,  It  Is  not  possible  to  render  more  than  a  specula¬ 
tive  account  of  the  meaning  of  our  measurements.  Bearing  this  In  mind, 
we  note  that  near  the  wall,  r/(D/2)  ■  0.75,  the  ratio  v'/u*  Increases 
persistently  between  X/D  ■  -1  and  e  ■  90®.  Because  this  coincides  with  a 
correspondingly  large  Increase  In  V,  we  suspect  that  the  term  -  vwaV/ar 
contributes  strongly  to  the  generation  of  v  (=  vc)  near  the  wall.  As 
of  e  ■  90°  the  ratio  v'/u'  decreases  near  the  wall.  The  decrease  Is  due 
to  a  gradual  reduction  In  v'  (as  opposed  to  an  Increase  In  u ' )  as  the 


flow  evolves;  see  Figs.  2  and  3.  Because  3V/3r  peaks  at  about  a  *  45°, 
from  this  position  forward  we  argue  that  -vw  3V/3r  now  contributes 
decreaslngly  to  the  balance  of  v'2.  This  then  allows  dissipation  and 
redistribution  mechanisms  to  reduce  the  level  of  v'2  near  the  wall*  as 
observed. 

The  core  flow,  r/D(2)  -  0,  experiences  a  marked  decrease  In  v'/u' 
between  a  ■  45*  and  90°  In  the  bend.  Between  these  locations  the  quan¬ 
tity  a V/30  is  large  and  positive.  Inspection  of  the  data  for  uv  obtained 
by  Chang,  et  al.  [5]  In  a  curved  duct  of  square  cross-section  suggests 
that  In  the  core  of  the  present  flowHy’  Is  positive.  As  a  result,  the 
term  -uv  aV/ae  should  be  expected  to  contribute  negatively  to  the  balance 
of  v  (and  hence  reduce  this  component)  between  e  »  45°  and  90°. 

The  Increase  In  v'/u'  between  e  »  90“  and  X/D  ■  1  In  the  core  of  the 
flow  Is  more  difficult  to  explain.  It  Is  probably  related  to  the  convec¬ 
tion  of  v'2  from  the  wall  periphery  (where  v'2  Is  generated)  into  the 
core. 

The  ratio  v'/u'  In  the  flow  entering  the  bend  at  X/D  ■  -1  Is  In 
qualitative  agreement  with  Laufer's  [193  pipe  flow  data,  both  at 
r/(D/2)  «  0  and  0.75.  In  the  downstream  tangent  the  expected  level  of 
v'/u'  is  re-established  near  the  pipe  wall  within  5  pipe  diameters.  By 
contrast,  the  turbulence  remains  strongly  anisotropic  In  the  core  of  the 
flow,.  The  persistence  of  this  anisotropy  may  be  related  to  the  presence 
of  an  extra  pair  of  symmetrical  vortical  structures,  discussed  earlier  in 
relation  to  cross-stream  flow  reversal. 
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3.3  Effect  of  Flow  Cross-Section 

A  comparison  was  made  between  some  of  the  present  results  for 
Re  •  57,400  and  corresponding  measurements  obtained  by  Chang,  et  al.  [5] 
In  a  180*  curved  duct  of  square  cross-section  with  Rc/D  ■  3.35  and 
Re  *  56,700.  Plots  of  square(s)-to-round(r)  ratios  of  circumferential 
mean  velocities  and  corresponding  turbulence  Intensities  are  shown  In 
Fig.  7. 

Absolute  values  of  CV/Ul)]s/CV/U|>]r  considerably  larger  than  unity 
Illustrate  clearly  the  dominant  Influence  exercised  by  the  flat  walls  on 
most  of  the  flow  through  a  curved  square  duct.  (Negative  values  in  the 
plots  are  due  to  the  cross-stream  flow  reversal  occurring  In  the  pipe  but 
not  In  the  duct.) 

At  X/D  ■  -1  the  Invlscld  transverse  motion  Induced  through  longitu¬ 
dinal  curvature  In  the  pipe  was  found  to  be  twice  as  strong  as  that  In 
the  duct.  However,  because  the  Re  and  Rc/D  ratios  between  the  two  con¬ 
figurations  Is  so  close,  potential  flow  theory  would  predict  essentially 
equal  cross-stream  flow.  Therefore,  we  conclude  that:  a)  the  more 
intense  secondary  motion  subsequently  Induced  in  the  duct  (due  to  the 
transverse  pressure  gradient)  acts  to  suppress  the  Initial  Invlscld  com¬ 
ponent  of  transverse  flow;  and,  b)  this  effect  roaches  further  upstream 
Into  the  straight  duct  than  Into  the  straight  pipe. 

The  circumferential  turbulence  Intensity  ratios  plotted  In  Fig.  7 
show  larger  fluctuation  levels  near  the  flat  wall  of  a  curved  duct  than 
near  the  side  wall  of  a  curved  pipe.  This  Is  In  keeping  with  the  higher 
velocities  and  hence  higher  rates  of  shear  expected  along  the  flat  walls 
of  a  curved  duct.  Surprisingly,  in  the  core  of  the  flow  the  fluctuations 
are  always  larger  In  the  pipe  and  may  be  related  to  the  additional  pair 


of  symmetrical  counter-rotating  vortices  which  develop  in  the  core  of  the 
curved  pipe.  The  relatively  large  value  of  [v'/U^lg/Cv'/Ujjilp  shown  at 
X/D  ■  -1,  near  the  wall,  Is  attributed  to  corner  contributions  to  fluc¬ 
tuations,  present  only  in  the  square  duct. 

4.  CONCLUSIONS 

This  study  places  on  a  firm  quantitative  basis  some  of  the  earlier 
qualitative  Interpretations  proposed  by  Rowe  [10]  for  the  mean  behavior 
of  turbulent  flow  through  a  180°  curved  pipe.  Specifically,  the 
existence  of  two  cross-stream  flow  reversals  Is  confirmed.  This  supports 
the  notion  of  an  additional  (symmetrical)  pair  of  counter-rotating  vor¬ 
tical  structures  embedded  In  the  core  of  the  flow  within  the  curved  pipe 
and  part  of  the  downstream  tangent.  Fully  developed  flow  conditions  are 
never  attained  In  the  present  bend. 

Interpretation  of  our  measurements  points  to  strong  generation  of 
circumferential  fluctuations  along  the  pipe  periphery  due  to  the  main 
secondary  flow.  By  contrast,  negative  contributions  to  generation  dampen 
circumferential  fluctuations  In  the  core.  It  Is  argued  that  the  dam- 
penlng  Is  offset  by  cross-stream  convection  of  v*  from  the  bend 
periphery  Into  the  core  of  the  flow. 

•A  comparison  between  the  cross-stream  motions  In  curved  circular 
pipes  and  curved  square  ducts  shows  that  the  latter  have  considerably 
higher  levels  of  secondary  flow.  Cross-stream  flow  reversals  were  not 
observed  In  the  square  duct  configuration.  It  is  suggested  that  the 
larger  relative  levels  of  cross-stream  fluctuations  detected  In  the  core 


of  a  curved  pipe  flow  are  probably  related  to  the  presence  there  of  an 
additional  pair  of  symmetrical  counter-rotating  vortices. 

Continued  work  should  focus  on  extending  the  data  base  for 
clarifying  the  role  of  the  additional  pair  of  vortical  structures  In 
curved  pipe  flow  and  the  relaxation  process  occurring  In  the  downstream 
tangent. 
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Figure  1.  The  test  section,  configuration  and  defini¬ 
tion  of  coordinate  system. 
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Figure  4.  Longitudinal  coiTiponents  of  mean  velocity 
and  turbulence  intensity  in  180°  curved 
pipe  flow  with  n.c/0  ■  3,375  and 
Re  °  110,000;  { - )  is  Laufer's  [19]  dat 
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5,  Circumferential  components  of  mean  velocity 
and  turbulence  Intensity  In  180°  curved 
pipe  flow  with  Rc/D  ■  3.375  and 
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APPENDIX  I 


Derivation  of  Beam  Location  and  Velocity  Correction  Formulae 

Nomenclature 

d  wall  thickness 
fjj  Doppler  frequency 

r  real  beam  Intersection  location  In  water  {relative  to  pipe  center) 
ry  virtual  beam  intersection  location  in  air  (relative  to  pipe  center) 

R  pipe  inner  radius 

U  "corrected"  longitudinal  (streamwlse)  velocity;  obtained  usings  -y  In 
U  ■  fD  x/(2  sin  i|i) 

UM  measured  longitudinal  velocity;  obtained  by  usings  ■  a  In 
UM  *  fD  x^2  s1n 

V  ''corrected"  circumferential  (cross-stream)  velocity;  obtained  by  using 
i|i  ■  y  In  V  ■  fD  x/(2  sin  if) 

VM  measured  circumferential  velocity;  obtained  by  usings  ■<*  In 

#  VM  ’  fD  X/(Z  s1n  ♦) 
a  half  angle  between  light  beams  In  air 

y  half  angle  between  light  beams  In  water 

refractive  Index  of  water  (relative  to  air) 

Up  refractive  Index  of  plexiglass  (relative  to  air) 

X  wavelength  of  laser  light 
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Longitudinal  Velocity  Component  -  Straight  and  Curved  Pipe  Sections 

(See  Figure  1-1) 


Assume  small  angles  throughout  such  that: 
sin  (angle)  «  tan  (angle)  *  angle 
From  Snell's  law, 
sin  a  ■  np  sin  p  . 

For  small  angles  a,  p  » 

3  “  «/np  .  (1) 

From  the  geometry  of  Fig*  1-1, 

a  ■  (R  +  d  -  ry)  tan  a  «  (R  +  d  -  ry)  a  , 
f  *  d  tan  p  ■  d  a/np  . 

From  Snell 's  law, 

sin  y  ■  sin  p  , 
nf 

With  (1)  and  a  small  angle  approximation  fort. 

Y  *  o/>.f  .  (2) 

Solving  for  r  from  geometry  along  with  the  above  relations  for  a,  f  andy? 

♦ 

ft  -  p  ■  (a  -f)/tan  y  *  a  [R  -  ry  +  d  (1  -  l/np)]  nf/a 
or, 

R  -  r  •  [R  -  ry  +  d  (1  -  l/np)3  nf  , 
p  >  CR  (l/nf  -  1)  +  ry  -  d  (1  -  l/np)3nf  , 


(3) 
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Because 

UM  ~  sTil  a  ' 

It  follows  that 

U  *  sin  a  UM/s1n  y  *  a  U^/y  . 

Then  from  (2): 

U“T1fuM  (4) 
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Fig.  1-1.  Effect  of  parallel  flat  surface-  on  beam  Intersection  location  and 
half  angle,,  not  to  scale. 


1-4. 


From  Snell's  law  and  small  angle  approximations: 

*3/np  , '  (5) 

<2  “  ®  rif /n p  •  (6) 

From  the  law  of  sines  and  small  angle  approximations: 

-  (R  +  d)/R  ,  (7) 

y/o  ■  R/r  .  (8) 


As  approximations  to  a,  with  small  angles: 

a  «  (R  +  d  -  ry)  o 
and 

a  ■  (R  +  d)  (a  “  $ )  • 

Equating  these  expressions  yields, 

3  -  ry  c/(R  +  d)  .  (9) 

Frpm  (5),  (6)  and  (7), 

o  ■  (R  +  d)  0/(R  flf)  .  .  (10) 

From  (9)  and  (10), 

o  »  ry  a/(R  nf)  .  (11) 


From  (8)  and  (11), 


nf  Y 

We  now  seek  an  expression  fory/a*  the  half  angle  correction  factor. 

From  geometrical  considerations,  we  may  use  the  following  approximations 
for  e: 


a  -  0  -  y  +  a 


'str nr  " R  » 


e  *•  d  6. 


Then,  from  the  above  two  expressions, 


o  -  p  -  y  +  o  *«-i 


Introducing  (5),  (9)  and  (11), 


rv“  .xrv“_dp  _d  rva 

a .  irrj  -  y  +  inrp  ir  ^  *  wrwrT~  • 

* 

Dividing  by  a, 

,  rv  .  ,  .  rv  _  d  rv 

1  “  rn  • Y  /a  +  Tnrjr  *  ■(R-^ajT'"'  * 


Solving  t'orY/os 
a  ■  1-n 


— -TirraT  +?t; 


or,  equivalently, 

Y_.rv  i  d 

a  nr  [  ~m-  m 1  ’  “jnf|r 


■]  +  1  , 


(12)  then  becomes 


Circumferential  Velocity  Component  -  Curved  Pipe  Section  (See  Figure  1-3) 


Snell's  law  and  small  angle  approximations  yield: 

&  -  «/np  ,  (15) 

♦  ■  <*>  np/nf  .  (16) 

From  the  geometry  of  Fig.  1-3  and  using  (15), 
a  ■  (R  +  d  -  ry)  a  , 
u  ■  5  “  o  » 

f  •  b  tan  0  ■  d  tan  0  ■  d  0  -  d  — 

nP 

g-a«f»[R+d(l-  l/np)  -  ry]  o  , 

■ 

and  with  the  preceding  expression  for  g, 

o  “  sin-1  g/R  »  a  [1  -  rv/R  +  d/R  (1  -  lAip)]  .  (17) 

From  (15)  and  (17)  substituted  Into  the  expression  forw, 

i 

«  -  O  [l/np  -  1  +  rv/R  -  d/R  (1  -  l/np)] 
or 

*<■»"«  [ry/R  -  (1  +  d/R)  (1  -  l/np)3  .  (18) 

Using  (16)  and  (18), 


hiii  '  Hi  ii  nfl 


h 


♦  ■fl[(i)]i*  (1  +  d/R)  ££-—)] 
w  n  f  n  {  n  { 


(19) 


From  the  geometry  of  Fig.  1-3, 

h  ■  r  sin  a  ■  (R  -  r  cos  a)  tan  4>  . 

Solving  for  r  with  small  angles  a,  lit 
r«R(jt/(t|>+o)  . 

This  becomes,  with  (17)  and  (19): 

R  [jj—  •  (1  +  d/R)  (^  -  l/nf)1 

r  ■  — JL— — ■  i I, . . . .  — - — — 


or 

-  R  Ejr  (1  +  d/R)  (np/nf  -  l/nf) 

+  1  +  d/R  (1  -  1/rip  +  1/n^  -  np/nf)  +  l/nf  -np^f 
/  C(1  -  np/nf)  i"  +  R  n p/n f ] 

Here  y  is  defined  by 

Y  «  o  +  ♦ 

and  can  be  found  using  (17)  and  (19). 

S^lnce  the  expression  for  V  is  given  by 

v  -v«  -TOT  *  VM y*  -VMTOT  • 


It  follows  that, 


nt 


R 


Fig.  I- 
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.  Effect  of  flat-to-curved  surfaces  on  beam  Intersection  location 
and  half  angle;  not  to  3ca1e. 


APPENDIX  II 


II. 1 


Frequency  Shifting  Procedure 

The  electronic  downmlxer  Is  one  of  two  modular  units  on  the  DISA  55 
Nll-16  frequency  shifting  device.  Either  one  of  two  channels  (A  or  B)  may  be 
used.  The  photomultiplier  (PM)  output  signal  cable  Is  attached  to  the  PM  IN 
port  In  the  rear  of  the  downmlxer  and  another  Is  attached  from  the  MIXER  OUT 
to  the  PM  SIGNAL  IN  port  of  the  DISA  counter.  The  least  complicated  mode  of 
operation  Is  to  set  the  net  electronic  shift,  f$,  greater  than  the  maximum 
frequency  expected  In  the  flow.  For  instance,  If  a  maximum  center  line  velo¬ 
city  of  5  m/s  Is  expected,  with  an  estimated  minimum  calibration  fringe  fcctor 
based  on  the  optical  parameters  and  test  section  geometry  of  10  m/s/MHz,  then 
a  simple  ratio  gives  a  0.5  MHz  upper  frequency  estimate. 

With  the  shift  frequency  denoted  by  fto  and  signal  frequency  by  fQ,  we 
would  then  set  the  shift  dial  corresponding  to  f|_Q  <  40  MHz  to  a  value 
slightly  greater  than  0.S  MHz,  say  700  KHz.  The  signal  from  the  counter, 
monltdred  on  the  oscilloscope,  now  has  a  net  positive  mixer  frequency  shift 
imposed.  In  the  computer  program  for  evaluating  velocity,  the  shift  value 
fLo  must  be  subtracted  from  the  frequency  measured  by  the  counter.  In  the 
present  example,  we  would  add  -0.7  MHz  to  each  measured  Doppler  count.  Values 
of  velocity  computed  In  this  manner  are  either  positive  or  negative, 
depending  on  the  direction  In  which  the  fringes  are  moving.  The  frlnge  move¬ 
ment  Is  set  by  turning  tha  Bragg  cell  screw. 

A  verification  of  fringe  sense  of  translation  Is  necessary,  A  rotating 
diffraction  grating,  with  the  laser  beams  intersecting  just  on  the  grating 
surface,  can  be  used  to  achieve  this.  Computed  values  of  velocity  and  their 


signs  for  prescribed  speeds  and  sense  of  grating  rotation  should  agree  with 
the  physically  (known)  Imposed  conditions.  Wherever  the  calculated  sign  of 
velocity  Is  positive  the  fringes  are  moving  opposite  to  the  grating  direc¬ 
tion;  the  opposite  Is  true  for  negative  values  of  velocity. 

If  a  mirror  Is  used  to  divert  the  beam  path,  as  In  the  present  case,  care 
must  be  taken  to  note  which  way  the  fringes  move  after  beam  reflection. 


Ill  .1 


APPENDIX  III 
DATA  TABULATION 


The  following  Is  a  listing  of  measurements  obtained  at  various  longitudi¬ 
nal  stations  in  the  test  section.  Scans  were  always  performed  along  the 
radial  direction.  In  the  plane*  «  w  /2  (see  Fig.  1).  Both  unweighted  and 
weighted  results  are  tabulated.  For  weighting  the  one-dimensional  scheme  pro¬ 
posed  by  McLaughlin  and  Tlederman  [17]  was  used.  Differences  between  weighted 
and  unweighted  data  are  discussed  In  the  text  where  It  Is  concluded  that 
weighting  the  circumferential  components,  V  and  v‘,  yields'  unrealistic 
results. 
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1.  19 

0.  050 

0.  027 

C.  041 

0.  000 

J 

1.  23 

0.  046 

C.  029 

0.  034 

1 

1.  20 

0.  046 

0.  026 

0.  039 

Weighted  Date 


■  37,400 


9  a  48 


Re  «  110,000 


r/tD/2>  I  U/Ub 


0.  967 

1 

0.  S3 

0.  123 

0.  224 

0.  119 

1 

0.  83 

0.  943. 

1 

0,  89 

0. Ill  . 

0.  193 

0.  11R 

1 

0.  89 

0.  077 

t 

0.  93 

0.  099 

0.  119 

0.  103 

1 

0.  94 

0.  010 

1 

0.  90 

0.  093 

0.  063 

0.  074 

1 

0.  96 

0.  742 

1 

0.  99 

0.  092 

0.  029 

0.  004 

1 

0.  97 

0.  673 

1 

1.  01 

0.  \:4}? 

0.  013 

0.  041 

1 

0.  99 

0.  607 

5 

1.  03 

0.  0S7 

0.  003 

0.  030 

1 

1.  01 

0.  340 

l 

1.  03 

0.  001 

-0.  003 

0.  036 

1 

1.  02 

0.  472 

I 

1.  07 

0.  001 

-0.  009 

0.  034 

1 

1.  04 

0.  403 

1 

1.  09 

,  0.  074 

0.  010 

0.  038 

1 

1.  06 

0.  337 

1 

1.  11 

0.  060 

-0.  023 

0.  040 

1 

1.  08 

0.  270 

l 

1.  12 

0.  063 

-0.  031 

0.  046 

1 

1.  09 

0.  202 

i 

1.  14 

0.  037 

-0.  037 

0.  047 

1 

1.  11 

0.  133 

1 

1.  16 

0.  033 

-0.  046 

0.  031 

1 

1.  12 

0.  067 

l 

1.  17 

0.  049 

-0.  061 

0.  034 

1 

1.  13 

0.  000 

1 

1.  10 

0.  047 

-0.  060 

0.  030 

1 

1.  14 

0.  122 

0.  226 

0.  098 

0.  109 

0.  204 

0.  107 

0.  009 

0.  123 

0.  093 

0.  004 

0.  033 

0.  068 

0.  080 

0.  032 

0.  035 

0.  081 

0.  017 

0.  041 

0.  078 

0.  007 

0.  032 

0.  073 

>0.  002 

0.  029 

0.  072 

-0.  007 

0.  030 

0.  070 

-0.  011 

0.  031 

0.  067 

-0. 016 

0.  024 

0.  062 

-0.  027 

0.  039 

0.  060 

-0. 032 

0.  042 

0.  036 

-0.  039 

0.  046 

0.  031 

-0. 033 

0.  032 

0.  040 

-0. 061 

0.  054 

Unweighted  Data 

Re  -  37,400 

Re  -  110,000 

r/(D/2>  I  U/Ub  u ' /Ob 


0.  967 
0.  943 
0.  077 


0.  742 
0.  673 


0.  340 


0.  405 


0.  270 
0.  202 
0.  133 
0.  067 
0.  000 


v'/Ub  I  U/Ub 


u'/Ub 


0.  121 

0.  286 

0.  096 

1 

0.  87 

0.  119 

0.  268 

0.  083 

0.  107 

0.  261 

0.  094 

1 

0.  90 

0.  103 

0.  257 

0.  083 

0.  093 

0.  199 

0.  113 

1 

0.  93 

0.  087 

0.  193 

0.  084 

0.  093 

0.  133 

0.  093 

1 

0.  96 

0.  083 

0.  130 

0.  080 

0.  091 

0.  080 

0.  083 

1 

0.  98 

0.  078 

0.  096 

0.  073 

0.  083 

0  ,  041 

0.  074 

1 

1.  00 

0.  078 

0.  063 

0.  074 

0.  086 

0.  017 

0.  073 

1 

1.  02 

0.  077 

0.  023 

0.  064 

0.  081 

-0.  014 

0.  066 

1 

1.  02 

0.  074 

-0. 003 

0.  038 

0.  080 

-0.  023 

0.  038 

1 

1.  04 

0.  071 

-0.  022 

0.  056 

0.  072 

-0.  049 

0.  058 

1 

1.  06 

0.  069 

-0.  037 

0.  052 

0.  067 

-0.  037 

0.  054 

1 

1.  08 

0.  066 

~0.  031 

0.  032 

0.  064 

-0.  074 

0.  065 

1 

1.  09 

0.  061 

-0.  066 

0.  049 

0.  057 

-0.  081 

0.  062 

1 

1.  11 

0.  039 

-0.  074 

0.  048 

0.  051 

-0.  093 

0.  033 

1 

1.  12 

0.  033 

-0.  083 

0.  031 

0.  048 

-0.  103 

0.  052 

1 

1.  13 

0.  031 

“0.  098 

0.  031 

0.  047 

-0.  112 

0.  037 

1 

1.  14 

0.  047 

~0.  106 

0.  031 

Weighted  Date 


Re  •  97# 400 


Re  *  HO,  000 


r/<D/2>  J  U/Ub 


U/Ub 


967 
0.  949 
877 
810 
0.  742 
0.  679 
0.  607 
0.  9*0 
0.  472 
0.  409 
0.  337 
0.  270 
0.  202 
0.  139 
0.  067 
0.  000 


0.  9 
1.  0 
1.  1 
1.  1 
1.  1 
1.  1 
1.  1 
1.  0 
1.  O 
1.02 
0.  98 
0.93 
0.  91 
0,  88 
0.  89 
0.89 


0.  167 

0.  099 

0.  104 

1 

0.  93 

0.  133 

0.  086 

0.  094 

? 

1.  03 

0.  092 

0.  043 

0.  069 

1 

1.  09 

0.  084 

0.  022 

0.  099 

1 

1.  09 

0.  070 

0.  011 

0.  046 

1 

* 

1.  08 

0.  000 

0.  002 

0.  042 

! 

0.  070 

-0.  005 

0.  041 

1 

0  C71 

-o.  009 

0,  040 

1 

0.  073 

-0.  009 

0.  039 

1 

0,  082 
0.  089 
0.  100 
0.  100 
O.  100 
0.  102 
0.  098 


-0, 010 
-0.  007 
-0.  005 
-0.  002 
-0.  OOC 
0.  001 
0.  002 


0.  040 
0.  037 
0.  037 
0.  036 
0.  036 
0.  039 
0.  034 


0.  90 
0.  94 
0.  90 
0.  86 
0.  83 
0.  31 
0.  82 


0.  129 

0.  083 

0.  086 

0.  109 

0.  069 

0.  079 

0.  078 

0.  031 

0.  057 

0.  074 

0.  016 

0.  043 

0.  066 

0.  008 

0.  036 

0.  066 

-0.  003 

0.  033 

0.  066 

-0.  003 

0.  032 

0.  066 

-o.  ooa 

0.  034 

0.  070 

-0.  011 

0.  034 

0.  077 
0.  093 
0.  110 
0.  116 
0.  119 
0.  121 
0.  123 


-0.  011 
-0.  009 
-0.  008 
-0.  003 
-0.  004 
“0.  003 
•*0.  002 


0.  033 
0.  033 
0.  032 
0.  031 
0.  032 
0.  031 
0.  031 


Unweighted 

Data 

Re  *  97, 

400 

Re  ^  110, 

000 

r/CD/2)  ! 

U/Ub  u'/Ub 

V/Ub 

v'/Ub 

t 

1 

u/Ub 

u  '/'Jil 

V/Ub 

v'/Ub 

0. 967  \ 

1. 

0 

2  0.  199 

0.  197 

0.  101 

1 

1.  00 

0.  119 

0.  163 

0.  083 

0. 945-  l 

1. 

0 

9  0.  123 

0.  176 

0.  0V6 

1 

1.  04 

0.  103 

0.  146 

0.  088 

0.  877  I 

t. 

1 

4  0.  089 

0.  113 

0.  094 

1 

1.  09 

0.  076 

0.  100 

0.  983 

0.810  *1 

1 

1 

4  0. 082 

0.  077 

0.  092 

{ 

1.  09 

0.  072 

0.  062 

0.  079 

0.  742  J 

1. 

1 

4  0.  076 

0  040 

0.  090 

V 

i 

l.  08 

0.  069 

0,  035 

0.  071 

0. 679  l 

1. 

1 

2  0.  073 

0.  010 

0.  083 

1 

t.  06 

0.  066 

0.  014 

0.  068 

0.  607  I 

1. 

1 

0  0.  069 

-0.  017 

0.  0G2 

1 

1.  09 

0.  066 

-0, 008 

0.  066 

0.  340  J 

t 

0 

8  0.  070 

-C.  029 

0.  073 

1 

S.  03 

0.  066 

-0.  027 

0.  069 

0.  4^2  1 

1. 

0 

5  0.  072 

-0.  036 

0.  074 

1 

01 

0.  070 

-0.  039 

0.  069 

0.  409  1 

1. 

0 

a  o  oai 

-0. 036 

0.  074 

1 

0.  98 

0.  076 

-0.  040 

0,  068 

0.  337  1 

0 

9 

8  0.  087 

-0.  028 

0.  071 

1 

0,  99 

0.  089 

-0,  036 

0.  064 

0.  270  I 

0, 

9 

3  0.  098 

-0.  024 

0.  072 

1 

0.  91 

0.  104 

-0.  032 

0.  067 

0.  202  1 

0 

9 

1  0.  098 

-0.  OJ 3 

0.  069 

1 

0.  87 

0.  112 

-0.  023 

0.  064 

0.  139  ! 

0 

8 

8  0.  098 

-0.  008 

0,  073 

1 

0.  83 

0.  118 

-0.  023 

0.  067 

0.  067  ! 

0 

8 

7  0.  102 

-0.  007 

0.  067 

f 

0.  S3 

0.  121 

-0.  017 

0.  064 

0,  000  1 

0. 

8 

7  0,  098 

-0.  003 

0.  066 

1 

0.  84 

0.  121 

-0.  015 

0.  066 

e  «  ia9° ' 


Weighted  Date 

Re  -  57 

,400 

Re  *  110 

,  000 

r/(D/2> 

t  U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

1 

U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

0.  967 

1  1.06 

0.  095 

0.  052 

0.  074 

1 

1.  06 

0.  086 

0.  044 

0.  074 

0.  945 

1  1.  09 

0.  084 

0.  040 

0.  074 

1 

1.  09 

0.  078 

0.  013 

0.  046 

0.  877 

1  1. 08 

0.  075 

0.  011 

0.  056 

1 

1.  09 

0.  068 

-0.  006 

0.  041 

0.  810 

1  1.  05 

0.  080 

-0.  012 

0.  054 

1 

1.  06 

0.  073 

-0.  017 

0.  047 

0.  742 

1  1.  02 

0,  082 

-0.  022 

0.  056 

1 

1.  04 

0.  077 

-0.  018 

0.  047 

0.  675 

1  1.01 

0.  082 

-0.  023 

0.  033 

1 

1.  01 

0.073 

-0.  014 

0.  043 

0.  607 

1  0.  99 

0.  081 

-0.  013 

0.  048 

1 

0.  99 

0.  077 

-0.  008 

0.  030 

0.  540 

1  0.  98 

0.  081 

-0.  009 

0.  042 

1 

0.  97 

O;  076 

-0.  008 

0.  033 

0.  472 

1  0.  97 

0.  080 

-0.  003 

0.  041 

1 

0.  96 

0.073 

0,  000 

0.  035 

0.  405 

I  0.  96 

0.081 

0.  004 

0.  039 

1 

0,  95 

0.  076 

0.  003 

0.  035 

0.  337 

1  0. 96 

0.  081 

0.012 

0.  041 

1 

0.  94 

0.  077 

o.  ooa 

0.  038 

0.  270 

1  0.  95 

0.  081 

0.  013 

0.  041 

1 

0,  94 

0.  079 

0.  010 

0.  039 

0.  202 

l  0.  95 

0.  083 

0.  013 

0.  042 

! 

0.93 

0.  083 

0.  008 

0.  040 

0.  135 

1  0.  95 

0.  082 

0.017 

0.  044 

1 

0.  93 

0.  080 

o.  ooa 

0.  039  * 

0.  067 

1  0.  95 

0.  084 

0.014 

0.  044 

1 

0.  93 

0.  083 

0.  007 

0.  039 

0.  000 

l  0.  93 

0.  084 

0.017 

0.  046 

1 

0.  93 

0.  079 

0.  006 

0.  041 

Unweighted  Date 

Re  -  57, 

400 

Re  »  110, 

000 

r/tD/2) 

1  U/Ub 

U'/Ub 

V/Ub 

v'/Ub 

1 

U/Ub 

u'/Ub 

V/Ub 

v.'/Ub 

0.  967 

<  1.07 

0.  092 

0.  127 

0.  092 

1 

1.  07 

0.  084 

0.  135 

0.  098 

0.  945 

1  1.  09 

0.  082 

0.  108 

0.  110 

1 

1.  09 

0.  077 

0.  032 

0.  089 

0. 877  # 

1  1.  09 

0.  074 

0.  037 

0.  110 

1 

1.  09 

0.  068 

-0.  030 

0.  083 

0.  810 

t  1.06 

0.  079 

-0.  041 

0.  101 

t 

1.  06 

0.  072 

-C. 039 

0.  OBO 

0.  742 

1  1. 03 

0.  031 

-0. 067 

0.  039 

! 

l.  04 

0,  076 

-0.  063 

0.  079 

0.  675 

1  1. 02 

0.  081 

“0.  068 

0.  082 

1 

1.  02 

0.  074 

™0.  033 

0.  080 

0.  607 

1  0.  98 

0.  081 

-0.  053 

0.  OBO 

1 

0.  99 

0.  076 

-0. 036 

0.  076 

0  540 

)  0. 98 

0.  079 

-0.  034 

0.  076 

1 

0.  98 

0,  075 

-0.  009 

0.  069 

p.  472 

1  0. 97 

0.  084 

-0.  019 

0,  072 

1 

0.  97 

0.  0/9 

-0.  006 

0.  073 

0.  405 

1  0.  96 

0.  081 

0.  009 

0.  071 

1 

0.  96 

0.  075 

0.  013 

0.  072 

0.  337 

J  0. 96 

0.  081 

0.  027 

0.  071 

1 

0.  95 

0.  077 

0.  023 

0.  072 

0.  270 

I  0. 96 

0.  081 

0.  032 

0.  068 

1 

0.  94 

0.  078 

0.  027 

0.  075 

0.  202 

1  0.  93 

0.  083 

0.  033 

0.  069 

1 

0,  94 

0.  082 

0.  026 

0.  072 

0.  135 

1  0.  98 

0.  082 

0,  033 

0,  070 

1 

0.  93 

0,  OSO 

0.  023 

0.  070 

0.  067 

1  0. 96 

0.  084 

0.  036 

0.  070 

1 

0.  94 

0,  083 

0.  020 

0.  075 

0.  000 

1  0.  95 

0.  081 

0.  037 

0.  074 

1 

0.  94 

0.  079 

0.  018 

0.  070 

I  ,  * 


0  -  177° 


W«lght«d  Data 


Ra  -  57 

o 

o 

** 

Ra  «  110 

,000 

r/(D/2> 

I  U/Ub 

U'/Ub 

V/Ub 

v'/Ub 

1 

U/Ub 

U'/Ub 

V/Ub 

v'/Ub 

0.  967 

1  0.  98 

0.  096 

0.  042 

0.  067 

i 

1.  02 

0.  081 

0.  032 

0.  060 

0.  949 

1  0,  99 

0.  081 

0.  019 

0.  093 

1 

1.  02 

0.  076 

0.  020 

0.  04? 

0.  877 

1  1.  00 

0.  067 

-C.  002 

0.  042 

f 

1.  02 

0.  068 

M7.  008 

0.  040 

0.810 

i  0.  99 

0.  061 

“0.  016 

0.  034 

f 

i.  oo 

0.  069 

-0.  016 

0.  037 

0.  742 

1  0. 98 

0.  099 

-o.  on 

0.  030 

1 

1.  00 

0.  061 

-0.  010 

0.  029 

0.  679 

i  0.98 

0.  099 

-0.  007 

0.  028 

1 

0.  99 

0.  060 

«0.  004 

0.  023 

0.  607 

1  0.  99 

0.  060 

-0.  004 

0.  028 

1 

0,  99 

0.  060 

-0,  001 

0.  029 

0.  940 

1  1.  00 

0.  062 

0.  001 

0.  028 

! 

1.  00 

0.  060 

0.  002 

0.  024 

0.  472 

1  1.  00 

0.  062 

0.  002 

0.  027 

i 

1.  00 

0.  060 

0.  004 

0.  023 

0.  409 

1  1,  00 

0.  067 

0.  003 

0.  029 

t 

1.  00 

0.  062 

0.  004 

0.  026 

0.  337 

1  1.  00 

0.  066 

0.  005 

0.  030 

i 

1.  00 

0.  063 

0.  005 

0.  026 

0.  270 

1  1.  01 

0.  069 

0.  009 

0.  030 

i 

0.  99 

0.  061 

0.  005 

0.  027 

0.  202 

!  1. 01 

0.  067 

0.  009 

0.  034 

i 

0.  99 

0.  063 

0.  004 

0.  029 

0.  139 

I  1.  OS 

0.  067 

0.  004 

0.  037 

j 

0.  99 

0.  069 

0.  004 

0.  030 

0.  067 

1  1.  01 

0.  069 

0.  007 

0.  039 

i 

0.  98 

0.  064 

0.  009 

0.  031 

0.  000 

t  1.  01 

0.  067 

0.  006 

0.  039 

f 

0.  98 

0.  063 

0.  006 

0.  032 

Unwaighbad  Data 

i 

Ra  «  37, 

400 

R*  *  110, 

000 

. .  — - 

r/(D/2> 

t  U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

i 

U/Ub 

uVUb 

V/Ub 

v'/Ub 

0.  967 

I  0  98 

0.  101 

0.  113 

0.  109 

i 

1.  02 

0.  080 

0.  106 

0.  087 

0,  949  „ 

1  0.99 

0.  081 

0.  096 

0.  093 

\ 

1.  02 

0.  070 

0.  076 

0.  086 

0.  877 

1  1. 00 

0.  067 

-0.  004 

0.  080 

i 

1.  02 

0,  070 

-0.  030 

0.  078 

0.010 

t  1.  00 

0.  061 

-0.  043 

0.  093 

i 

1.  00 

0.  064 

-0.  098 

0  058 

0.  742 

3  0.  99 

0.  060 

-0.  032 

0.  048 

i 

1.  00 

0.  099 

-0.  039 

0.004 

0.  679 

1  0.99 

0.059 

-0.  022 

0.  047 

i 

0.  99 

0.  099 

-0,  01? 

0.  050 

0.  607 

\  0.  99 

0.  060 

•*0.  Oil 

0.  049 

i 

0.  99 

0.  058 

-0.  009 

0.  032 

0.  940 

:  i.  oo 

0.  063 

-0.  003 

0.  048 

i 

1.  00 

0.  062 

“0  002 

0.  048 

6.  472 

1  1.  01 

0.  062 

0.  000 

0.  049 

i 

1.  00 

0.  061 

0.  006 

0.  030 

0.  409 

t  1.  00 

0.  067 

0.  004 

0.  091 

i 

1.  00 

0.  061 

0.  008 

0.032 

0.  3.J7 

1  1.  01 

0.  066 

0.  010 

0.  049 

1.  00 

0.  0«3 

0.  010 

0.  091 

0.  270 

1  1. 01 

0.  069 

0.  009 

0.  061 

1 

0.  99 

0.  061 

o.  oil 

0.  033 

0.  202 

J  1. 01 

0.  067 

0.  010 

0.  068 

1 

0.  99 

0.  063 

0.  008 

0.  057 

0.  135 

1  1. 02 

0.  067 

0.  006 

0.  076 

1 

0.  98 

0  066 

0.  006 

0.  062 

0.  067 

i  1.  01 

0.  066 

0.  011 

0.  06S 

1 

0.  99 

0.  064 

0.  011 

0.  062 

0.  000 

I  1.  01 

0.  067 

0.  012 

0.  067 

1 

0.  99 

0.  062 

0.  012 

0.  060 

Weighted  Date 

■ 

R#  -  57,400 

Re  -  110.000 

r  ■■  "  *-  1 

r/<D/2>  1  U/Ub 

u'/Ub  V/Ub 

v'/Ub 

1  U/Ub  •  u'/Ub  V/Ub 

v'/Ub 

0.  967 

1 

0.  94 

0.  105 

0.  028 

0.  046 

1 

0.  96 

0.  138 

0.  026 

0.  049 

0.  949 

1 

0.  98 

0.  100 

0.  024 

0.  048 

1 

0.  98 

0.  084 

0.  018 

0.  041 

0.  877 

t 

1.  02 

0.  104 

0.  006 

0.  032 

1 

0.  99 

0.  082 

0.  004 

0.  027 

0.  810 

t 

1.  01 

0.  084 

0.  001 

0.  028 

1 

1.  00 

0.  061 

-0.  010 

0.  030 

0.  742 

1 

1.  03 

0.  067 

-0. 002 

0.  026 

i 

1.  00 

0.  056 

-0.  01 1 

0.  028 

0.  675 

l 

1.  02 

0.  074 

-0.  001 

0.  026 

1 

0.  98 

0.  056 

-0.  003 

0.  021 

0.  607 

1 

1. 03 

0.  067 

0.  002 

0.  025 

1 

0.  99 

0.  094 

-0.  003 

0.  021 

0.  540 

1 

1.  04 

0.  058 

0.  007 

0.  027 

1 

1.  00 

0.  056 

0.  000 

0.  022 

0.  472 

i 

1.  04 

0.  059 

0.  008 

0.  029 

1 

1.  00 

0.  057 

0.  004 

0.  022 

0.  405 

i 

1.  05 

0.  062 

0,  007 

0.  031 

l 

0.  99 

0.  059 

0.  005 

0.  023 

0.  337 

i 

1.  05 

0.  062 

0,  007 

0.  030 

1 

0.  99 

0.  061 

0.  005 

0.  026 

0.  270 

i 

1.  04 

0.  062 

0.  008 

0.  030 

1 

0.  99 

0.  061 

0.  007 

0.  027 

0.  202 

i 

1.  04 

0,  062 

0.  008 

0.  031 

1 

1.  00 

0.062 

0.  008 

0.  029 

0.  135 

i 

1.  05 

0,  064 

0,  009 

0,  031 

1 

1.  00 

0.  063 

0.  007 

0.  028 

0.  067 

i 

1.  05 

0.  064 

1 

1.  00 

0.  062 

0.  009 

0.  030 

0.  000 

i 

1.  05 

0.  064 

1 

1.  00 

0.  062 

0.  007 

0.  030 

Unweighted  Date 

Re  *  57.400 

R*  m  110.000 

r/(D/2>  I  U/Ub 


0.  967 
0.  945 
0.  877 
0.  810 
0.  742 
0.  675 
0.  607 
0.  540 
0.  472 
0.  405 
0/337 
0.  270 
O'.  202 
0.  135 
0.  067 
0.  000 


U/Ub 


0.  98 

0. 

113 

0.  072 

0.  072 

0.  99 

0. 

130 

0.  074 

0.  077 

1.  02 

0. 

119 

0.  020 

0.  058 

1.  02 

0. 

094 

0.  002 

0.  052 

1* 

1.  02 

0. 

080 

-0.  008 

0.  047 

1.  02 

0. 

093 

-0.  007 

0.  060 

1.  03 

0. 

081 

0.  004 

0.  049 

1.  04 

0. 

098 

0.  016 

0.  053 

1.  05 

0. 

059 

0.  018 

0.  094 

1.  05 

0. 

061 

0.  014 

0.  063 

1.  09 

0. 

062 

0.  017 

0.  057 

1.  05 

0. 

061 

0.  020 

0.  050 

1.  09 

0. 

062 

0.  019 

0.  054 

1,  05 

0. 

064 

0.  021 

0.  058 

1.  09 

0. 

064 

1.  05 

0. 

063 

0.  190 

0.  081 

0.  075 

0,  091 

0.  065 

0.  070 

0.  090 

0.  017 

0.  056 

0.  061 

-0. 033 

0.  091 

0.  059 

-0. 039 

0.  045 

0.  055 

-0. 013 

0.  042 

0.  055 

-0. 014 

0.  044 

0.  055 

-0.  003 

0.  041 

0.  056 

0.  011 

0.  042 

0.  058 

0.  014 

0.  043 

0.  060 

0.  019 

0.  091 

0.  060 

0.  021 

0.  050 

0.  061 

0.  023 

0.  051 

0.  062 

0.  021 

0.  053 

0.  062 

0.  026 

0.  054 

0,  062 

0.  029 

0.  056 

x/d  ■  +2 


Weighted  Data 

Ra  ■  37* 

400 

• 

Rt  -  110*000 

r/<D/2> 

1 

U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

l 

U/Ub 

U'/Ub 

V/Ub 

v'/Ub 

0.  967 

0.  06 

0.  132 

0.  020 

0.  040 

1 

0.  93 

0.  103 

0.014 

0.  037 

0.  943 

0.  92 

0.  090 

0.  010 

0.  036 

1 

0.  93 

0.002 

0.010 

0.  033 

0.  077 

1.01 

0.  066 

0.  010 

0.  024 

1 

1.  02 

0.  063 

0.  006 

0.  023 

0.  010 

1.  03 

0.  032 

0.  003 

0.  023 

1.  01 

0.037 

-0.  001 

0.  021 

0.  742 

1.  02 

0.032 

0.  004 

0.  024 

1 

0.  99 

0.  032 

-0.  002 

0.  020 

0.  673 

1.02 

0.  049 

0.  003 

0,  022 

1 

0.  99 

0.032 

-0.  002 

0.  021 

0.  607 

1.  02 

0.  030 

0.  003 

-  0.  024 

1 

1.  00 

0.  033 

0.  001 

0.  021 

0.  340 

1.  02 

0.  033 

0.  003 

0.  020 

1 

1.  01 

0.  036 

0.  003 

0.  022 

0.  472 

1.  02 

0.  036 

0.  007 

0.  029 

» 

1.  01 

0.  030 

0.  003 

0.  023 

0.  403 

1.  03 

0.  030 

0.  000 

0.  029 

i 

1.  01 

0.060 

0.  000 

0.  023 

0.  337 

1.  03 

0.  060 

0.  009 

0.  031 

1 

1.  01 

0.  060 

0.  00B 

0.  027 

0.270 

1.03 

0.  061 

0.  010 

0.  032 

1 

1.01 

0.  062 

0.009 

0.  030 

0.  202 

1.  04 

0.  062 

0.  000 

0.  030 

1 

1.  02 

0.  064 

0.  on 

0.  022 

0.  133 

1.  04 

0.  063 

0.  010 

0.  032 

1 

1.  01 

0.  062 

0.  010 

0.  031 

0.  067 

1.  03 

0.  063 

0.  014 

0.  034 

1 

1.  01 

0.  062 

0.  010 

0.  031 

0.  000 

1.  03 

0.  062 

0.  012 

0.  032 

1 

1.  01 

0.  062 

0.  007 

0.  029 

Unwttghtad  Data 


Ra  -  87* 

400 

Rt  ■  110* 

000 

r/(D/2) 

I  U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

1 

U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

0.967 

1  0.  00 

0.  122 

0.  036 

0.  063 

1 

0.94 

0.  090 

0.  033 

0.  064 

0.  943 

1  0. 93 

0.  093 

0.  04a 

0.  036 

1 

0.  90 

0.  073 

0.  040 

0.  063 

0.  077  1  , 

1  1.  01 

0.  066 

0.  026 

0.  037 

1 

1.  02 

0.  063 

0.  023 

0.  030 

0.  010 

1  1.  03 

0.  032 

0.  016 

0.  043 

1 

1.01 

0,  037 

-0.  004 

0.  041 

0.  742 

1  1.  02 

0.  033 

0.  010 

0.  041 

! 

1.  00 

0.  033 

-0.  007 

0.039 

0.673 

1  1.  02 

0.  049 

0.  006 

0.  039 

I 

0.  99 

0.  052 

-0.  000 

0.  040 

0.  607 

1  1.02 

0.  030 

0.  003 

0.  041 

1 

1.  00 

0.  033 

0.  000 

0.  041 

0.  340 

.1  1.02 

0,  033 

0.  012 

0.  032 

1 

1.  01 

0.  033 

0.  003 

0.  041 

0.  472 

1  1. 02 

0.  037 

0.  010 

0.  031 

1 

1.  02 

0.  030 

0.  on 

0.  041 

0.  403 

1  1.  04 

0.  030 

0.  020 

0.  032 

1 

1.  03 

0.  060 

0.  024 

0.  046 

0.  337 

i  1.04 

0.  060 

0.  021 

0.  067 

1 

1.01 

0.  060 

0.  024 

0.  030 

0.  270 

1  1.04 

0.  061 

0.  024 

0.  063 

1 

1.  02 

0.  062 

0.  030 

0.  032 

0.  202 

1  1.04 

0.  062 

0.  019 

0.  061 

1 

1.  02 

0.  064 

0.  033 

0.  030 

0.  133 

1  1. 04 

0.  064 

0.  026 

0.  034 

1 

1.  02 

0.  062 

0.  032 

0.  054 

0.  067 

1  1. 03 

0.  063 

0.  034 

0.  056 

1 

1.  01 

0.  062 

0.  032 

0.  034 

0.  000 

1  1.  03 

0.  061 

0.  033 

0.  032 

1 

1.  02 

0.  062 

0.  026 

0.  034 

x/D  -  +3 


41 

4 


■] 

l 

;j  Weighted  Data 


p 

Re  -  97 

,  400 

Re  -  110 

,  000 

1 

r/<D/2> 

1  U/Ub 

u  VUb 

V/Ub 

vVUb 

1 

U/Ub 

u  VUb 

V/Ub 

‘.7 

0.  967 

1  0.  93 

0.  076 

0.  019 

0.  039 

1 

1.  00 

0.  073 

0.  002 

0.  031 

ki 

'1 

0.  949 

1  0.  99 

0.  069 

0.  017 

0.  039 

1 

1.  02 

0.  067 

0.  009 

0.  029 

0.  877 

1  0.  98 

0.  093 

0.  on 

0.  030 

1 

1.  02 

0.  097 

0.  002 

0.  023 

1 

0.  810 

1  1.  00 

0.  048 

0.  009 

0.  026 

1 

1.  01 

0.  089 

-0.  001 

0.  023 

0.  742 

1  1.  00 

0.  046 

0.  009 

0.  024 

1 

1.  02 

0.  093 

-0.  001 

0.  021 

0.  679 

1  1.  02 

0.  092 

0.  009 

0.  024 

1 

1.  03 

0.  097 

0.  000 

0.  021 

0.  607 

1  1.  02 

0.  092 

0.  003 

0.  024 

1 

.... 

0.  002 

0.  022 

.. 

0.  940 

1  1.  02 

0.  099 

0.  004 

0.  029 

! 

1.  03 

0.  098 

0.  004 

0.  0213 

, 

0.  472 

1  1.  02 

0.  097 

0.  003 

0.  027 

1 

— 

0.  007 

0.  026 

i 

0.  409 

t  1.  02 

0.  061 

0.  004 

0.  027 

1 

1.  04 

0.  060 

0.  009 

0.  028 

0.  337 

1  1.  02 

0.  069 

0.  006 

0,  029 

1 

1.  04 

0.  069 

0.  011 

0.  031 

0.  270 

1  1.  02 

0.068 

0.  007 

0.  030 

1 

1.  04 

0.  068 

0.  012 

0.  031 

• 

0.202 

1  1.  03 

0.  067 

0.  003 

0.  031 

1 

1.  04 

0.  068 

0.  013 

0.  033 

0.  139 

1  1.  02 

0.  068 

0.  008 

0.  031 

l 

1.  03 

0.  068 

0.  013 

0.  033 

ii 

0.  067 

t  1.  02 

0.  068 

0.  007 

0.  031 

1 

1.  02 

0.  068 

0.  010 

0.  032 

i 

0.000 

1  1.  02 

0.  067 

0.  003 

0.  031 

l 

1.  02 

0.  071 

0.  009 

0.  031 

►  ‘’1 

Unweighted  Data 

i' 

L' 

Re  -  97, 

400 

Re  -  110, 

000 

i 

r/(D/2) 

1  U/Ub 

u  VUb 

V/Ub 

v  VUb 

1 

U/Ub 

u  VUb 

V/Ub 

v  VUb 

0.  967 

1  0.  94 

0.  074 

0.  049 

0.  064 

1 

1.  01 

0.  072 

0.  030 

0.  063  ' 

\ 

0.  949 

i  0.  99 

0.  067 

0.  048 

0.  061 

1 

1.  02 

0.  069 

0,  019 

0.  092 

ij 

0.  877 

I  0.  99 

0.  093 

0.  033 

0.  049 

1 

1.  02 

0.  098 

0.  002 

0.  048 

0.  810 

J  1.  00 

0.  048 

0.  019 

0.  049 

1 

1.  02 

0.  096 

-0. 006 

0.  048 

i 

0.  742 

1  1.  00 

0.  047 

0.  013 

0.  041 

1 

1.  02 

0.  093 

-0.  004 

0.  041 

m 

0.  679 

\*  1.01 

0.  049 

0.  014 

0.  040 

l 

t.  04 

0.  097 

0.  000 

0.  044 

:* 

0.  607 

i.  1.02 

0.  092 

0.  006 

0.  041 

1 

— 

0.  007 

0.  043 

■ 

0.  940 

>  1.  02 

0.  099 

0.  011 

0.  044 

l 

1.  04 

0.  098 

0.  014 

0.  047 

V 

0.  472 

1  1.  02 

0.  097 

0.  008 

0.  047 

1 

— 

0.  023 

0.  048 

0.  409 

1  1.  02 

0.  061 

0.  011 

0.  048 

5 

1.  04 

0.  060 

0.  030 

0.  092 

| 

0.  337 

I  1.  02 

0.  064 

0.  014 

0.  093 

\ 

1.  04 

0.  064 

0.  037 

0.  054 

r 

0/270 

>  1.  02 

0.  067 

0.  018 

0.  092 

1 

1.  04 

0.  063 

0.  039 

0.  099 

0.  202 

1  1.  03 

0.  066 

0.  022 

0.  093 

« 

1.  04 

0.  068 

0.  042 

0.  054 

0.139 

1  1.  03 

0.  067 

0.  022 

0.  092 

1 

1.  04 

0.  068 

0.  042 

0.  059 

!>; 

0.  067 

I  1.  03 

0.  066 

0.  020 

0.  094 

1 

1.  03 

0.  063 

0.  037 

0.  060 

0.  000 

1  1.  03 

0.  067 

0.  022 

0.  032 

1 

1.  03 

0.  070 

0.  029 

0.  060 

X/D  ■  +4 


If! 


',lA 

>'V. 

r/3 


I 


Uaightcd  Data 


R«  «  97*  400 


R*  *  1 10* 000 


—  — — — — —  — — —  — — — — —  —  <4i 


-;]  r/CD/2> 

1  U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

1 

U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

2  0.  967 

1  0.93 

0.  109 

0.  010 

0.  037 

1 

0.  94 

0.  084 

0.  004 

0.  029 

I  0.  949 

1  0.  93 

0.  081 

0.  010 

0.  036 

1 

0.  98 

0.  078 

0.  006 

0.  031 

f,  0.  877 

1  1.02 

0.  071 

0  006 

0.  027 

} 

1.  01 

0.  062 

0.  003 

0.  023 

R  0.  810 

1  1.09 

0.  064 

0.  003 

0.  024 

1 

1.  01 

0.  060 

0.  000 

0.  020 

0.  742 

1  1. 04 

0.  067 

0.  002 

0.  022 

1 

1.  01 

0.  093 

0.  000 

0.  020 

h  0.679 

t  1.  03 

0.  097 

“0.  001 

0.  021 

1 

1.  02 

0.  094 

-0.  001 

0.  020 

L  0.  607 

(  1.  09 

0.  097 

0.  C02 

0.  022 

1 

1.  01 

0.  062 

-0.  008 

0.  020 

I  0.  940 

i  1.  09 

0.  061 

0.  002 

0.  023 

l 

1.  02 

0.  062 

0.  002 

0.  022 

‘  0.  472 

t  1.  06 

0.  063 

0.  003 

0.  029 

1 

1.  04 

0.  063 

0.  004 

0.  023 

5  0.  409 

1  1.  07 

0.  072 

0.  009 

0.  027 

1 

1.  09 

0.  067 

0.  OOC 

0.  024 

;  0.  337 

1  1.  07 

0.  071 

0.  009 

0.  023 

1 

1.  04 

0.  069 

0.  006 

0.  029 

■I  0.  270 

1  1. 07 

0.  079 

0.  007 

0.  029 

1 

1.  03 

0.  071 

0.  006 

0.  023 

\  0.202 

t  1.  07 

0.  073 

0.  007 

0.  030 

l 

t.  04 

0.  070 

0.  008 

0.  027 

0,  139 

«  1.  06 

0.  080 

0.  007 

0.  030 

l 

1.  04 

0.  073 

o.  ooa 

0.  030 

0.  067 

t  1.  06 

0.  080 

0.  010 

0.  031 

l 

1.  04 

0.  072 

0.  009 

0.  030 

0.  000 

1  1.  06 

0.  081 

0.  010 

0.  031 

1 

1.  04 

0.  073 

0.  009 

0.  032 

Unw*ight*d  Dat« 

R*  ■  97, 

400 

R>  -  110* 

000 

r/(D/2> 

t  U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

1 

U/Ub 

u'/Ub 

V/Ub 

v'/Ub 

0.  967 

1  0.  99 

0.  IIS 

0.  029 

0.  069 

1 

0.  99 

0.  033 

0.  021 

0.  063 

0.  949 

l  0.  93 

0.  097 

0.  031 

G.  066 

1 

0.  98 

0.  079 

0.  029 

0,  064 

0.877 

L  t.  04 

0.  076 

0.  017 

0.  049 

l 

1.  01 

0.  063 

0.  013 

0.  049 

0.810 

!  1.  09 

0.  083 

0.  009 

0.  042 

1 

1.  01 

0.  062 

0.  003 

0.  040 

0.  742 

}  1.  09 

0.  081 

0.  006 

0.  043 

\ 

1.  01 

0.  094 

0.  000 

0.  038 

0.  679 

t  1.  04 

0.  064 

-0.  001 

0.  037 

l 

1.  02 
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FOREWORD 


This  technical  report  provides  a  summary  of  the  accomplishments 
relating  to  Contract  Number  N00013-80-C-QQ31  for  the  project  entitled  * 
‘‘Turbulent  Flow.and  Heat  Transfer  In  Passage  Around  130°  Bend  -  An 
Experimental  and  Numerical  Study."  The  project  was  funded  by  the  Office 
of  Naval  Research  under  the  supervision  of  Mr.  K.  Elllngsworth, 

i  ■  ( 

The  work  at  the  University  of  California,  Berkeley,  was  performed 
In  close  collaboration  with  Professor  Brian  E.  Launder  and  his  research 
group  at  the  University  of  Manchester  Institute  of  Science  and  Technology. 
Work  at  UCB  emphasized  measurements  and  computations  in  a  180n  bend  bf 
square  cross-section.'  Limited  measurements  were  also  obtained  1h  a  180° 
bend  of  circular  cross-section.  These  new  results  are  Included  in  the 
present  report.  Work  at  liMIST  has  emphasized  measurements  and  computations 
In  a  180°  bend  of  circular  cross-section.  In  addition,  measurements  were 
also  obtained  In  a  duct  of  square  cross-section.  The  work  at  OMlST  Is 
reported  elsewhere  by  Professor  Launder. 

While  the  studies  at  UCB  and  at.  UMIST  progressed  along  separate 
and  distinct  routes,  efforts  were  continuously  focussed  on  performing 
complementary  research  activities  and  sharing  knowledge  at  all  levels 
as  soon  as  It  was  available.  This  Interaction  hay  strengthened  and 
stimulated  the  respective  research  proqrains.  As  a  result.,  considerable 
pregrass  has  been  made  towards  measuring,  understanding  and  making 
predictable  the  flow  through  strongly  curved  passages.  However,  as  this 
report  will  explain,  much  work  remains  to  be  done. 
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PROJECT  SUMMARY 


A  brief  summary  Is  provided  below  of  the  most  significant  achievements 

relating  to  the  work  performed' for  this  contract.  The  summary  Is  divided 

into  two  parts.  .Section  LI  evaluates  the  work  performed  in  a  curved 

passage  of  square  cross-section,  while  Section  1.2  evaluates  more  recent 

Information  obtained  In  a  curved  passage  of  circular  cross-section.  Where 

appropriate,  reference  Is  made  to  publications  already  available. 

/• 

1,1  Work  In  a  180°  Curved  Passage  of  Square  Cross-Section 

This  work  has  been  reported  In  complete  detail  In  the  Final  Annual 
Report  for  Period  October  1,  1982  -  September  30,  1983  (see  reference  2 
In  Section  2.2  of  this  report).  A  brief  summary  of  the  main  accomplishments, 
findings  and  recommendations  listed  In  that  report  Is  given  here. 

Main  Findings 

1.  Measurements  of  mean  velocities  and  their  associated  turbulent 
stresses  along  the  streamwlse  and  radial  directions  were  made  at 
several  longitudinal  stations  using  the  laser-Doppler  veloclmetry 
technique.  One  shear  stress  component  was  also  measured. 

2..  The  mean  flow  data  revealed  features  which  were  In  qualitative 
agreement  with  results  expected  from  Invlscld  flow  analysis. 

.3.  The  turbulent  stresses  displayed  previously  undocumented 
anisotropic  characteristics  arising  from  shearing  motions  induced 
In  the  core  of  the  flow. 

4.  Measurements  In  the  downstream  .tangent  showed  drastic 
reductions  In  the  bend-induced  secondary  motion  taking  place  In 
less  than  five  hydraulic  diameters. 


5.  Numerical  predictions  of  the  flow  using  a  semi-elliptic 
procedure  and  a  ie-e  model  of  turbulence  yielded  results  in  poor 
agreement  with  the  measurements.  False  physical  diffusion  and 
the  Isotropic  characteristics  Inherent  In  the  model  were  thought 
to  be  the  causes  for  the  differences  observed. 

6.  'Predictions  wlth'an  algebraic  stress  model  also  yielded 
results  In  poor  agreement  with  the  measurements.  However,  numerical 
experiments  showed  that  the  algebraic  stress  model  Is  capable  of 
resolving  the  anisotropy  of  the  turbulence  provided  that  the  mean 
velocity  field  Is  known  accurately. 

Among  the  recommendations  evolving  from  the  study  of  flow  through 
strongly  curved  passages  of  square  cross-section  were  the  following. 
Recommendations 

1.  More  detailed  measurements  are  required  to  understand  the 
process  of  flow  relaxation  taking  place  In  the  tangent  downstream 
of  the  bend. 

2.  Grid  refinement  techniques  and/or  higher  order  schemes  should 
be  explored  to  Improve  special  resolution  of  numerical  predictions 
of  the  flow  through  strongly  curved  passages. 

3.  While  the  use  of  an  algebraic  stress  model  is  recommended  for 
predicting  anisotropic  turbulent  flow  through  bends,  the  success  of 
this  closure  approach  relies  on  a  good  spaclal  resolution  of  the 
mean  flow. 

4.  It  should  be  established  what  are  appropriate  law-of-the-wall 
relations  for  turbulent  flows  through  strongly  curved  passages  of 
square  cross-section. 


